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Abstract 



We give a review and present a comprehensive calculation for the leptonic 
constant fs^ of the low-lying pseudoscalar and vector states of Bc-meson in 
the framework of static and QCD-motivated nonrelativistic potential models 
taking into account the one-loop and two-loop QCD corrections in the short 
distance coefficient that governs the leptonic constant of Be quarkonium sys- 
tem. Further, we use the scaling relation to predict the leptonic constant of 
the nS-states of the be system. Our results are compared with other models 
to gauge the reliability of the predictions and point out differences. 



I. INTRODUCTION 

The Be mesons provide a unique window into the heavy quark dynamics. Although they 
are intermediate to the charmonium and bottomonium systems the properties of Be mesons 
are a special case in quarkonium spectroscopy as they are the only quarkonia consisting of 
heavy quarks with different flavours. Because they carry flavour they cannot annihilate into 
gluons so are more stable and excited Sg-states lying below BD (and BD* or B*D) threshold 
can only undergo radiative or hadronic transitions to the ground state Be which then decays 
weakly. There are two sets of S-wave states, as many as two P-wave multiplets (the IP and 
some or all of the 2P) and one D-wave multiplet lying below BD threshold for emission of B 
and D mesons. As well, the F-wave multiplet is sufficiently close to threshold that they may 
also be relatively narrow due to angular momentum barrier suppression of the Zweig allowed 
strong decays. However, the spectrum and properties of these states have been calculated 
various times in the framework of heavy quarkonium theory [1-18]. 

The discovery of the Be meson by the Collider Detector at Fermilab (CDF) Collaboration 
[19-23] in the channel Be — J/iplu {I — e^ii^r) with low-lying level of pseudoscalar mass 
Mb, = QA^oaI GeV and hfetime tb, = 0.46lgi^ ± 0.03 ps confirmed the theoretical predic- 
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tions regarding various Be meson properties, spectroscopy, production and decay channels 
[1-18]. 

The description of various aspects of Be meson physics [24-27] has received recently a lot 
of attention both from theoretical as well as experimental sides. It could offer a unique probe 
to check the perturbative QCD predictions more precisely and lead to a new information 
about the confinement scale inside hadrons. The higher 6c-meson state cascades down 
through lower energy be state via hadronic or electromagnetic transitions to the pseudoscalar 
ground state Be meson. The theoretical study of the pure leptonic decays of the Be meson 
such as Be — > li^i can be used to determine the leptonic decay constant /s^ [28-33] , as well as 
the fundamental parameters in the standard Model (SM). Hence, it is the only meson within 
SM which is composed of two nonrelativistic heavy quarks of different flavors: open charm 
and bottom quarks. Its spectroscopy production mechanisms and decays differ significantly 
from those of charmonium J/%1) and upsilon T as well as hadrons with one heavy quark. 
Indeed, this meson is also a long lived system decaying through electroweak interactions. It 
stands among the families of cc and bb and thus could be used to study both quantitatively 
and conceptually existing effective low energy frameworks for the description of bound state 
heavy quark dynamics, hke NRQCD [34,35], pNRQCD [36-38] and vNRQCD [39,40]. 

The calculation of the Be leptonic decay constant Jb^ can be carried out either using 
QCD-based methods, such as lattice QCD [41], QCD sum rules [42-45], or adopting some 
constituent quark models [1-18]. So far, lattice QCD has only been employed to calculate 
the Be purely leptonic width. As the QCD sum rules, the Be leptonic constant, as well 
as the matrix elements relevant for the semileptonic decays were computed. Moreover, the 
leptonic constant has also been calculated by using the nonrelativistic (NR) potential model 
[1-18]. 

Ikhdair et al. [46-47] applied the NR form of the statistical model to calculate the 
spectroscopy, decay constant and some other properties of the heavy mesons, including the 
be system, using a class of three static quarkonium potentials. Davies et al. [41] predicted 
the leptonic constant of the lowest state of the Be system with the recent lattice calculations. 
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Eichten and Quigg [2] gave a more comprehensive account of the decays of the Be system 
that was based on the QCD- motivated potential of Buchmiiller and Tye [48] . Gershtein et 
al. [1] also published a detailed account of the energies and decays of the Be system using 
a QCD sum-rule calculations [42-45]. Pulcher [3,4] also calculated the leptonic constant of 
the Be quarkonium system using the treatment of the spin-dependent potentials to the full 
one-loop level and thus included effects of the running coupling constant in these potentials. 
Furthermore, he also used the renormalization scheme developed by Gupta and Radford 
[49-52]. Ebert et al. [13] have also calculated the leptonic constant of the pseudoscalar and 
vector Be quarkonium system using a relativistic model and then compared their calculation 
with the NR model. Kiselev et al. [53-57] calculated the leptonic constant of the pseudoscalar 
meson Be in the framework of QCD-motivated potential models taking into account the one- 
loop [53], the two-loop [54] correction matching and scaling relation (SR) [55-57]. Capstick 
and Godfrey [58] predicted the result of /b^ = 410 MeV for the decay constant of the Be 
meson using the Mock-Meson approach or other relativistic quark models. Ikhdair and Sever 
[59-64] have calculated the decay constant in the nonrelativistic and semi-relativistic quark 
model using the shifted large-N expansion technique. Finally, Godfrey [65] has calculated 
the spectroscopy of Be meson in the relativixzed quark model. Motyka and Zalewski [66] 
proposed a nonrelativistic Hamiltonian with plausible spin dependent corrections for the 
quarkonia below their respective strong decay threshold. They found the decay constants 
for the ground state and for the first excited S-state of he quarkonium to be Jb^ — 435 MeV 
and Jb^ = 315 MeV, respectively. 

In this study, we have considered the quarkonium potential models, which are usually 
used in the literature [1-18,59-64], in calculating the Be leptonic constant based on the NR 
potential model, (cf. e.g. [59-64] and references contained there). Further, we have applied 
these potential models, with a great success, to fit the entire heavy quarkonum spectroscopy 
[59-64,67-71]. We insist upon strict flavor- independence of its parameters. We also use 
the potential models to give a simultaneous account of the properties of the cc, hh and he 
quarkonium systems. 
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The contents of this paper are as follows: In section II, we present briefly the solution 
of the Schrodinger equation for the be quarkonium mass spectrum using the shifted large- 
N-expansion technique (SLNET). In section III, we introduce the necessary expressions for 
the one-loop and two-loop corrections to the leptonic constant available at present. Further, 
in section IV we present the phenomenological and the QCD-motivated potentials used in 
the present work. Finally, section V contains our calculations and remarks made for the 
leptonic constant in our approach and scaling relation. 



II. WAVE EQUATION 

In this section we consider the N-dimensional space Schrodinger equation for any spherically 
symmetric potential V{r). If ■0(r) denotes the Schrodinger 's wave function, a separation of 
variables 'ijj{r) = Yi^rn{(^, (p)u{r)/r^^'^^/^ gives the following radial equation (in units h — c — 
1) [59-64,67-73] 



[ 4/x dr^ 16//r^ j 

where = (^91^92) /(^gi + ^92) the reduced mass for the two interacting particles. 
Here En,i denotes the Schrodinger binding energy and k — N + 21 — a, with a representing a 
proper shift to be calculated later on and I is the angular quantum number. We follow the 
shifted or 1/k expansion method [67-73] by defining 

and also the energy eigenvalue expansion [59-71] 

m=0 ^ 

—1/2 

where x — k (r/ro — 1) with vq being an arbitrary point where the Taylor expansion is 

2 

being performed about and Q is a scale to be set equal to k at the end of our calculations. 
Inserting equations (2) and (3) into equation (1) gives 



1 I 1 (k (2-a) ^ (l-a)(3-a) ^| 



4A; 



m=0 



^2 oo 



(2/) = UrXnriy), 



(4) 



where the final analytic expression for the 1/k expansion of the energy eigenvalues appro- 
priate to the Schrodinger particle is [59-64,67-71] 



^ m=0 



(5) 



where rir is the radial quantum number. Hence, we formulate the SLNET (expansion as 
for the nonrelativistic motion of spinless particle bound in spherically symmetric potential 
V(r). The resulting eigenvalue of the N-dimensional Schrodinger equation is written as 
[59-64,67-71] 



^Ur — k 



+ 



1 ^riV{ro) 



16// Q 
1 \{l-a){3-a) 
k 16// 



+ 



1, (2-a) 



+ a 



(1) 



+ 



a 



(2) 



k 



2 ' 



(6) 



where a^^^ and a*^^-* are two useful expressions given by Imbo et al [72,73]. Comparing 
equation (5) with (6) yields 

Q 



Eo = V(ro) + 



16/ir^ 



2 ' 



(7) 
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27 8/x 



(8) 



Eo — ^ 
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Q 



(9) 
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(10) 



The quantity Tq is chosen so as to minimize the leading term. Eg [59-64,67-71], that is. 
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dro 

Therefore, Tq satisfies the relation 



and 



drl 



> 0. 



(11) 



Q = 8/iroV(ro), 



(12) 



and to solve for the shifting parameter a, the next contribution to the energy eigenvalue 
El is chosen to vanish [72,73] so that the second- and third-order corrections are very small 
compared with the leading term contribution. The energy states are calculated by consid- 
ering the leading term Eq and the second-order and third-order corrections, it implies the 
shifting parameter 



a = 2 - {2nr + 1) 



3 + 



1/2 



(13) 



Therefore, the Schrodinger binding energy to the third order is 



^l-a)(3-a) 
16a( 



a 



(2) 



k 



O 



(14) 



Once the problem is collapsed to its actual dimension (A^ = 3), one is left with the task 
of relating the coefficients of our equation to the one-dimensional anharmonic oscillator in 
order to read the energy spectrum. For the N — 3 physical space, the equation (1) restores 
its three-dimensional form. Thus, with the choice k = \fQ which rescales the potential, we 
derive an analytic expression that satisfies Tq in equations (12)-(13) as 



1 + 2/ + {2nr + 1) 



3 + 



ror^(ro) 
^'(ro) 



1/2 



8K^'(ro) 



1/2 



(15) 



where = 0, 1, 2, ■ ■ ■ stands for the radial quantum number and / = 0, 1, 2, • ■ ■ stands 
for the angular quantum number. Once Tq is being determined through equation (15), the 
Schrodinger binding energy of the g^g2 system in equation (14) becomes relatively simple and 
straightforward. We finally obtain the total Schrodinger mass binding energy for spinless 
particles as 



^(?i?2) = ra-q^ + 2E^^i. 
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(16) 



As stated before in [59-73] , for any fixed n the computed energies become more convergent as 
I increases. This is expected since the expansion parameter 1/N or 1/k becomes smaller as 
/ becomes larger since the parameter k is proportional to n and appears in the denominator 
in higher-order correction. 

On the other hand, the spin dependent correction to the nonrelativistic Hamiltonian, 
which is responsible for the hyperfine splitting of the IS-state mass level is generally used 
in the form (cf. e.g. Refs. [2,59-64,74-76])^ 

A£;HP = ^|^|^,s(o)r. (17) 

Like most authors (cf. e.g. [2,59-64,74-76]), we determine the coupling constant as{mc) 
from the well measured experimental hyperfine splitting of the IS(cc) state value 

AEnF = Mj/^ - M^, = 1 17 ± 2 MeV. ( 18) 

We use the value of the coupling constant to reproduce the spin-averaged data (SAD) or 
center-of-gravity (c.o.g.) of the lowest charmonium state value, that is, M^(IS). In order to 
apply this formula one needs the value of the wave function at the origin, this is obtained 
by solving the Schrodinger equation with the nonrelativistic Hamiltonian and the coupling 
constant Q;s(/x). The radial wave function at the origin [74-76] is determined by 

|fl.(0)r = 2^(^^). (19) 

where |i?is(0)|^ = 47r |'j/'is(0)|^ . Hence, the total mass of the low-lying pseudoscalar i?c-state 
is given by the expression 

Mb, (0-) = me + mb + 2Ei,o - 3A£;hf/4, (20) 



^At present, the only measured splitting of raS'-levels is that of r]c and J/ip, which allows us to 
evaluate the so-called SAD using M^(IS) = {SMj/^ + Mr]c)/^ and also M(nS) = My(nS) - 

{Mj/^-Mri,)/An [53,74-76]. 
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and also for the vector 5*-state 



Mb.(1-) 



(21) 



Finally, the square-mass difference can be simply found by using the expression 



M|,(l-) - M|^(0-) = 2AEuF K + mb + lE^fl - A£;hf/4] . 



(22) 



III. LEPTONIC DECAY CONSTANT OF THE Bc-MESON 



The study of the heavy quarkonium system has played a vital role in the development of the 
QCD. Some of the earliest applications of perturbative QCD were calculations of the decay 
rates of charmonium [77-80]. These calculations were based on the assumption that, in the 
NR limit, the decay rate factors into a short-distance (SD) perturbative part associated with 
the annihilation of the heavy quark and antiquark and a long-distance (LD) part associated 
with the quarkonium wavefunction. Calculations of the annihilation decay rates of heavy 
quarkonium have recently been placed on a solid theoretical foundation by Bodwin et al. 
[81]. Their approach is based on NRQCD, and effective field theory that is equivalent to 
QCD to any given order in the relative velocity v of the heavy quark and antiquark [82]. 
Using NRQCD [34,35] to seperate the SD and LD effects, Bodwin et al. [34] derived a general 
factorization formula for the inclusive annihilation decay rates of heavy quarkonium. The 
SD factors in the factorization formula can be calculated using pQCD, and the LD factors 
are defined rigorously in terms of the matrix elements of NRQCD that can be estimated 
using lattice calculations. It applies equally well to S-wave, P-wave, and higher orbital- 
angular-momentum states, and it can be used to incorporate relativistic corrections to the 
decay rates. 

In the NRQCD [83] approximation for the heavy quarks, the calculation of the leptonic 
decay constant for the heavy quarkonium with the two-loop accuracy requires the matching 
of NRQCD currents with corresponding fuU-QCD axial-vector currents [54] 
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= -Ai>cV^ and = &7^75C, (23) 

NRQCD -^"^'^ QCD ' ' V ^ 

where 6 and c are the relativistic bottom and charm fields, respectively, x| and ■^c are the 
NR spinors of anti-bottom and charm and is the four- velocity of heavy quarkonium. The 
NRQCD [83] lagrangian describing the i^c-meson bound state dynamics is 

>Cnrqcd = Aight + {iDo + DV(2me)) i^c + xl (^^o - DV(2mb)) + • • • , (24) 

where Cnght is the relativistic lagrangian for gluons and light quarks. The two-component 
spinor field ipc annihilates charm quarks, while Xb creates bottom anti-quarks. The relative 
velocity v of heavy quarks inside the ^c-meson provides a small parameter that can be used 
as a nonperturbative expansion parameter. To express the decay constant fs^ in terms of 
NRQCD matix elements we express in terms of NRQCD fields -^c and Xb- Only 

QCD 

the A = current-component contributes to the matrix element in the rest frame of the 
i?c-nieson is: 

{0\h'75cmP)) = ifB.P\ (25) 

where \Bc{P)) is the state of the meson with four-momentum P. It has the standard 
covariant normalization 

^ / rsM^BMd'p = 2E6^'\p' - p), (26) 

and its phase has been chosen so that Jb^ is real and positive. In NRQCD, the matching 
yields 

67%5C = Koxlipc + K2iBxby.r>iPc + (27) 

where Kg — Ko{mc,mb) and K2 — K2{mc,mb) are Wilson SD coefficients. They can be 
determined by matching perturbative calculations of the matrix element (0| 67^750 \ B^ that 
is mainly resulting from the operator Xjyi^c in 



67%5C |5c) = Ko (0| \Bc) + (0| (Dx.)^D^, \B,) 



QCD 



NRQCD 
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+ (28) 

NRQCD 



where the matrix element on the left side of (28) is taken between the vacuum and the state 
\B^ . Hence, equation (28) can be estimated as 

2 3Mb, 



TT 



\Rism 



(29) 



Onishchenko and Veretin [83] calculated the matrix elements on both sides of equation (28) 
up to a1 order. In one-loop calculation, the SD-coefficients are 

1 



Kn = 1 and Kq = 



"8At2' 



(30) 



with 11 defined after equation (1). Further, Braaten and Fleming in their work [84] calculated 
the perturbation correction to Kq up to order- ag (one-loop correction) as 

tts(Ai) 



i^o = 1 + Ci- 



TT 



(31) 



with ci being calculated in Ref. [84] as 



ci = - 



mb-rric . rrih 
2 In — 



TTib + rric rric. 

Finally, the leptonic decay constant for the one-loop calculations is 



(32) 



(l-loop) 
Be 



TT 



2 m — 

nib + TTT'c f^c 



J Br- 1 



(33) 



where the NR leptonic constant [85,86] is 

f NR _ 
JBe — 



TlMf 



l^is(0)| 



(34) 



and II is any scale of order or rric of the running coupling constant. On the other hand, the 
calculations of two-loop correction in the case of vector current and equal quark masses was 
made in [87,88]. Further, Onishchenko and Veretin [83] extended the work of Ref. [87,88] 
into the non-equal mass case. They found an expression for the two-loop QCD corrections 
to Sc-meson leptonic constant given by 



+ 



TT 



TT 



(35) 
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where C2(M///) is the two-loop matching coefficient and with ci,2 are explicitly given in Eq. 
(32) and (Ref. [83]; see equations (16)-(20) therein), respectively. In the case of Sp-meson 
and pole quark masses (m^ = 4.8 GeV, rric = 1.65 GeV), they found 



(.(2— loop) 



1 - 1.4J 



- 24.24 



TT 



asifUb 



TT 



NR 
Be ■ 



(36) 



Therefore, the two-loop corrections are large and constitute nearly 100% of one-loop correc- 
tion as stated in Ref. [83]. 



IV. SOME POTENTIAL MODELS 

The potential models are suitable for the phenomenological studies, because they can re- 
produce the model formulae or numbers for the quantities, used as input values (the level 
masses, for example). Hence the potential models can be considered as phenomenological 
meaningful fittings of some experimental values, but they can not restore a true potential, 
that does not exist due to the nonperturbative effects [55-57]. 



A. Static potentials 

It is easy to see that most phenomenological static potentials used in the Scrodinger 
equation in [59-64] may be gathered up in a general form [59-64,89] 

V{r)^-ar-'' + br'^ + Vo 0<a,(3<l, a,b>0, (37) 

where Vq may be of either sign. Here we assume that the effective ^^52 potential consists of 
two terms, one of which, VV(r) = — ar~", transforms like a time-component of a Lorentz 4- 
vector and the other, Vs{r) = br^, like a Lorentz scalar. These static quarkonium potentials 
are monotone nondecreasing, and concave functions which satisfy the condition 

V'{r) > and V"{r) < 0. (38) 

At least ten potentials of this general form, but with various values of the parame- 
ters, have been proposed in the literature (cf. [59-64] and references contained there). 
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The generality (37) comprises the following five types of potentials used in the litera- 
ture: Some of these potentials have a = /3 as in the Cornell {a — (3 = 1); Song- Lin 
(a = /3 = I) and Turin (a = /3 = |) potentials with same sets of fitting parameters 
used in our previous works [59-64]. Further, Song [90] has also used a potential with 

q; = /3 = |. 

On the other hand, potentials with a ^ 13 have also been popular as 



1. Martin potential 

The phenomenological power-law potential proposed by Martin (cf. e.g. [55-57,59-64]) 
has q; = 0, /3 = 0.1 of the form 

Vmir) = buiAMrf-' + cm, (39) 
is labeled as Martin's potential [55-57] with a given set of adjustable parameters 

[bu, Cm, Am] = [6.898 GeV^ \ -8.093 GeV, 1 GeV] , (40) 

and quark masses 

[me, rribl = [1.800 GeV, 5.174 GeV] . (41) 
(potential units are also in GeV). 

2. Logarithmic potential 

A Martin's power-law potential turns into the logarithmic potential of Quigg and Rosner 
[55-57,59-64] corresponds to a — P ^ and it takes the general form 

VL(r) = 6Lln(ALr) + CL, (42) 

with an adjustable set of parameters 

[bh, cl, Al] = [0.733 GeV, -0.6631 GeV, 1 GeV] , (43) 
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and quark masses 

K, rribl = [1.500 GeV, 4.905 GeV] . (44) 

The potential forms in (39), and (42) were also used by Kiselev in [55-57]. Further, they 
were also been used for ip and T data probing 0.1 fm < r < 1 fm region [53]. 

Further, Motyka and Zalewski [89] used a nonnrelativistic potential with a = 1 and 
/3 = I for the bb quarkonium and then applied it to the cc and be quarkonium systems. 
Grant, Rosner, and Rynes [91] have suggested a = 0.045, /3 = 0. Heikkila, Tornquist and 
Ono [92] tried a = 1, /? = |. Some very successful potentials known from the literature are 
not of this type. Examples are the Indiana potential [93] and the Richardson potential [94]. 

B. QCD-motivated potentials 

We use two types of the QCD-motivated potentials: The Igi-Ono (10) and an improved 
Chen-Kuang (CK) potential models. The details of these potentials can be traced in our 
previous works in [59-64]. 

V. NUMERICAL RESULTS AND CONCLUSIONS 

Based on our previous works [17], we determine the position of the charmonium center- 
of-gravity (c.o.g.) M^(IS) mass spectrum and its hyperfine splittings by fixing the coupling 
constant as{mc) in (17) for each central potential. Further, we calculate the corresponding 
low-lying (c.o.g.) Mx(lS) and consequently the low-lying Mb^(IS).^ In Table I, we estimate 
the radial wave function of the low-lying state of the be system, so that 

|i?is(0)| = 1.18 - 1.24 GeV^/^ (45) 



^Kiselev et al. [15] have taken into account that AMt(IS) = ^^AM0(1S) with as(T)/as('0) ^ 
3/4. On the other hand, Motyka and Zalewski [29] also found - ^l/l^- 
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for the set of the central potentials given in section IVA. Further, in Table II we present 
our results for the NR leptonic constant f^f = 466 ± 19 MeV and = 464 ± 19 MeV 
as an estimation of the potential models without the matching. Our results are compared 
with those of Gershtein, Likhoded and Slabospitsky [95-97], who used Martin's potential 
and with those of Jones and Woloshyn [98] . Moreover, the one-loop correction /^^"'""^^ and 
the two-loop correction fs'^"""^^ are also given in Table III. Therefore, in the view of our 
results, our prediction for the one-loop calculation is 

and also for two-loop calculation 

Our estimate for f^f' is fairly in good agreement with the estimates in the framework of 
the lattice QCD result [8], /^f = 440 ± 20 MeV, QCD sum rules [42-45], potential models 
[1-18,59-64] and also the scahng relation [55-57]. It indicates that the one-loop matching in 
Ref. [84] provides the magnitude of correction of nearly 12%. However, the recent calculation 
in the heavy quark potential in the static limit of QCD [53] with the one-loop matching [54] 
is 

Therefore, in contrast to the discussion given in [54], we see that the difference is not crucially 
large in our estimation to one-loop value in the Sc-meson. 
Our final result for the two-loop calculations is 

/g-i-p) = 315^1^ MeV, (49) 

the larger error value in (49) is due to the strongest running coupling constant in Cornell 
potential. 

Slightly different additive constants is permitted in this sector for a charmed mesons 
to bring up data to its (c.o.g.) value. However, with no additive constant to the Cornell 
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potential [99] , we notice that the smaller mass values for the composing quarks of the meson 
leads to a rise in the values of the potential parameters which in turn produces a notable 
lower value for the leptonic constant as seen in Tables II and III. 

On the other hand, for the QCD motivated two types 10 potential, our calculations 
for the ground state radial wave functions are presented in Table IV. In Table V, the NR 
leptonic constant of the pseudoscalar Be state are in the range 354 — 426 MeV, and for the 
Bl are 353 - 424 MeV, for the type I. Further, they are found 353 - 457 MeV and 351 - 454 
MeV, for the type II. For instance, we may choose f^f = 391 MeV with = 200 MeV, 
for the type I, and = 396 MeV with Ajjg = 300 MeV, for the type II. These results 
surpass the matching procedure as well as methods followed in other references, cf. Refs. 
[41,53,54]. Further, for the CK potential, we present the decay constant in Table VI. 

The scaling relation (SR) for the S-wave heavy quarkonia has the form [55-57] 

fn_ ( Mr,{bc) y f mc + mA ^ d 
M^ibc) \M,ibc)) \ Apt J n' ^ ^ 

where rric and are the masses of heavy quarks composing the i?c-meson, /i is the reduced 
mass of quarks, and o? is a constant independent of both the quark flavors and the level 
number n. The value of d is determined by the splitting between the 2S and IS levels or 
the average kinetic energy of heavy quarks, which is independent of the quark flavors and 
n with the accuracy accepted. The accuracy depends on the heavy quark masses and it is 
discussed in [55-57] in detail. The parameter value in (50), d ~ 55 MeV, can be extracted 
from the experimentally known leptonic constants of ijj and T [55-57]. Hence, in the view 
of Table VII, the SR gives for the IS-level 

f§^) ~ 4441^3 MeV, (51) 

for all static potentials used. Kiselev estimated /^^ — 400 ± 45 MeV [54] and f^^^ — 
385 ± 25 MeV [55-57], Narison and Chabab found f^^^ = 400 ± 25 MeV [100,101]. 

Overmore, we give the leptonic constants for the excited nS-levels of the be quarkonium 

/or) \ 

system in Table VII. We remark that the calculated value of Jb^^Is) ~ MeV is in 
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good agreement with /b^(2S) — 280 ± 50 MeV being calculated in [53] for the 2S-level in the 
be system and it is also consistent with the scaling relation [55-57]. 

Finally, we have also noted that the leptonic constant is practically independent of the 
total spin of quarks, so that 

fv,n ^ fp,n = fn, (52) 

where Mb^ ~ m;, + rric- Thus, one can conclude that for the heavy quarkonium, the QCD 
sum rule approximation gives the identity of leptonic constant values for the pseudoscalar 
and vector states. 

In this work, we have successfuly applied the SLNET using a class of static and QCD- 
motivated potentials to calculate numerically the leptonic constant of the pseudoscalar and 
vector Bc-Taeson. Once the experimental leptonic constant of the Bf, meson becomes clear, 
one can sharpen the analysis. 
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TABLES 



TABLE I. The characteristics of the wave function |V'is(0)|^ and the radial wave function 
\Ris{0)\ = 47r|V'is(0)r both at the origin (in GeV^) obtained from the Schrodinger equation for 
various potentials. 



Level 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


CornelP 




0.112 


0.123 


0.111 


0.119 


0.102 


0.065 


\Ris{o)f 


1.413 


1.54 


1.397 


1.495 


1.28 


0.814 



^Here we cite Ref. [98] for the fitted set of parameters. 



TABLE IL Pseudoscalar and vector decay constants (/p = /b^, fv = fs*) of the Be meson, 
calculated in the different potential models (the accuracy is 5%), in MeV. 



Constants 


Cornell^ 


Song-Lin 


Turin 


Martin 


Logarithmic 


Cornell'' 


[95] 


[97] 


fNR 


464.5 


485.1 


462.0 


478.1 


441.7 


351.5 


460±60 


420±13 


fNR 


461.5 


482.2 


459.2 


475.3 


439.3 


349.7 


460±60 





'^Here Vq ^ 0. 

''HereVb = 0, see Ref. [98]. 



TABLE III. One-loop and two-loop corrections to pseudoscalar and vector decay constants of 
the low-lying he system, calculated in the different potential models (the accuracy is 4%), in MeV. 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


Cornell 


j.{\—loop) 
J Be 


393.6^ 


424.4 


399.6 


421.2 


399.3 


311.9 


j,{2—loop) 
J Be 


264.1^ 


333.0 


296.6 


339.2 


340.9 


238.1 


n{l — loop) 

Jb* 


391.0 


421.9 


397.1 


418.7 


397.1 


310.3 


f{2—loop) 

Jb* 


262.4 


331.0 


294.8 


337.2 


339.0 


236.9 



*To the first correction for all potentials, Kq = 0.85 — 0.90. 
''To the second correction for all potentials, Kq = 0.57 — 0.77. 
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TABLE IV. The characteristics of the radial wave function at the origin |i?is(0)|^ (in GeV^) 
obtained from the Schrodinger equation for the Igi-Ono potential model. 



LeveP 


lUU 


zOU 


dOO 


400 


500 


Type I*^ 














0.066 


0.08 


0.092 


0.095 


0.089 


\Rism' 


0.826 


1.005 


1.156 


1.19 


1.114 


Type IF 














0.065 


0.071 


0.082 


0.096 


0.109 


\Rism' 


0.819 


0.891 


1.03 


1.204 


1.366 



^For the fitted parameters, see Ref. [59]. 

•"We shifted the cc spectra by c = -22 to -31 MeV. 

'^We shifted the cc spectra by c = —15 to —26 MeV. 
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TABLE V. Pseudoscalar and vector decay constants of the Be meson, calculated in the Igi-Ono 

polcuilial model, in McV. 

















inn 


9nn 


^nn 


4nn 


5nn 


lype i .as — 


n 1 QQc 
u.±yy 


n 01 7 


n 0^8 

U.Zoo 


u.zo 






354.1 


391.1 


420.0 


426.0 


411.7 


fNR 

Jb* 


352.6 


389.2 


417.7 


423.6 


409.4 


Type W^-.as = 


0.199 


0.227 


0.230 


0.241 


0.251 


fNR 
JBc 


352.7 


368.2 


396.1 


428.6 


456.9 


fNR 

Jb* 


351.2 


366.4 


394.1 


426.4 


454.4 



^For the fitted parameters we cite Ref. [59] . 
'^We shifted the cc spectrum by c = —22 to —31 MeV. 
'^Input value AEhf= 117 MeV; determines as{mc)- 
^We shifted the cc spectra by c = —15 to —26 MeV. 



TABLE VL The radial wave function, pseudoscalar and vector decay constants of the Be meson, 
calculated in the improved Chen-Kuang potential model, in MeV. 



Quantity^ 


100 


180 


300 


350 


375 


400 


450 


500 


1^15(0)1^ 


1.017 


1.017 


1.017 


1.017 


1.017 


1.015 


0.957 


0.829 


fNR 
JBc 


393.5 


393.5 


393.5 


393.4 


392.5 


389.9 


381.8 


353.8 


fNR 

Jb* 


391.4 


391.4 


391.4 


391.3 


390.4 


387.8 


379.8 


352.3 



^Input value AEhf= 117 MeV; determines as(mc) = 0.27. 
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TABLE VII. The nS-levels leptonic constant of the be system, calculated in the different po- 
tential models (the accuracy is 3 — 7%), in MeV, using the SR. 



Quantity 


Cornell 


Song-Lin 


Turin 


Martin 


Logarithmic 


/l5 


449.6 


450.4 


448.0 


448.8 


420.9 


/2s 


305.8 


305.0 


303.3 


303.5 


284.7 


/35 


243.0 


243.2 


241.3 


241.8 


227.2 


f4S 


206.0 


207.1 


204.9 


205.9 


193.8 
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